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CORRIGENDA. 

Page 12. In the second line in the proof of HerscherB theorem, insert 

.-. /(e«) = &c. 

Page 16. The Ex. in Art. 32, will not be correct, except 6 and not 1 
be taken for the increment. 

Page 60. The complementary function may be further simplified into 

** ^^KB <■ • • • • » 



2 DIFFERENCES. 

troduced this variable, and performed the operation, regarding 
unity for its increment, we can express the result in terms 

of X by writing 7 for z. 



Differences of Explicit Functions. 

8. To find the difference of Au^-¥ B; u, being a function 
of Xy and A and B quantities which do not depend on x. 

A(Au^ + JB) = (-4w,,i + JB) - (Au^ + £) = -4w,„ - Au^ = AAu^, 

Hence when -4 « 0, we have AB = 0. 

4. To find the difference of a number of functions of x, 
connected by the operations of addition or subtraction. 

= Au^ + At?„ - Aw. 

m g g 

5. To find the difference of the product of two functions. 

6. To find the difference of the quotient of two functions. 

A ?^ = ^' - ^ = "»>■". ~ ^'«"« 

_ («, + A«,) », - (t>, + At),) M, 
».Am, - M At). 



s m 
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DIFFERENCES. 

7. To find the difference of the continued product of any 
number of successive values of a function. 

Let uu^^ .... u^^^ be the product. 

Hence if u^^ A^-*- A^x, 

8. To find the difference of the quotient of continued pro^ 
ducts of any number of successive values of two functions u^ 
and «?,. 

Let * '*^ .' . ' " — Si? , be the quotient. 



V^i .... t?,,^ t?,^^<),^ ^>«»<^ ^'A.i ^. 



«+i» 






t^-.i • • . • ti. 



Hence where the numerator is unity, and Vg^^A^ + Ax^ 

9. To find the differences of any rational integral function. 
Let aa^ + baf''^ + . . • . Aa? + / »= w, be the function. Then 

»r waa;""' + b^af^ + , . . . t,a: + A^, 

1^2 



4 DIFFBRBNCBS. 

which is a rational integral function lower by one dimension 
than that from which it is derived. 

Repeating the operation^ and denoting the result or A (At^,) 
by A'ti,* ^^ htLve 

A'ti, »».(»- 1) aaj""' + bji^ + . ^ . . hjc ^i^^ 

and so on^ and finally^ when we have repeated the operation (n) 
times^ using the same notation^ 

A^tsfj, « n .(n - 1) . . . . d.2.1a. 

Hence the n^ difference is constant^ and therefore all 
differences of higher orders vanish. 



10, To find the differences of sin u^ and cos u,. 
A sin ti, « sin u,^^ - sin w, = sin (u, + Aw J - sin w. 



. . Aw 
2 8m — 



Aw, / Aw,\ 

cos w, + — -' I , 

2 \ ' 2 )' 



= 2 sin — ^ sin {u, + J(ir + Awj}, 

A cos tl, = cos tl^^j - cos U, = COS {Ux + AwJ - COS u 



= - 2 sm 



Aw, . / AwA 
— = sin Ug + — - I , 

= 2 sin — ' cos (w, + J(ir + Aw,)}. 
Making w, « aa; + /3, and repeating the operation (n) times, 
A"sin(cui; + /3) = f 2sin^ j sin<(eui:+ /3) + - (ir -f a)l, 

A" cos (aic + j3) = ( 2 sin ^ J cos {(fix + j3) + - (ir + a)}. 



mFFERENCES. 

11. To find the differences of tan w,, and tan"'w,. 
A tan u, = tan (w, + AwJ - tan «,, 

sin (Ug + Aw,) cos u^ - cos (w, + Aw J sin w. 



cos w^^j cos u. 



sin Aw, 



cos w,^j cos w. 



Whence, A tan ax = 



cos (a: + 1) acosoo;' 
A tan'^w, = tan*^w,^j - tan'^w,, 

^M.i -^s _ ._-! Aw, 



= tan-' -^ ^ = tan 



1 + w,,,w. 1 + w,,,w. 



• 



Whence, A tan"* ax = tan .^ — ; -r — x . 

•1 + (a? + l)a;a* 

12. To find the differences of a", 

Ao-^a'^'-a'^Ca- l)a*. 
Repeating the operation, 

AV = (a - lja% 
and so on, till we come to 

AX = (« - !/«•. 



13. To find the difference of log w,. 

Alog w, = log w,^j - logw, = log 



w 









b BEPARATION OF SYMBOIi^. 

Sueeessive Vdluei of a Function* 

14. If in any function of x, as «,, we substitute successively 
the values « + A, x ^ 2h . .. . x -¥ nh, for x, the results are 
termed the successive values of the function. 

The first successive value is thus expressed^ Du, « u^^j^ ; and 
h is generally unity, except where the contrary is stated. 

Hence, I/u,'- Du^^^=^ u,^^, ^w, - «,*3A> &c., 

and generally, i>X * w#*«a • 



iMcUions between the Successive Vdl/ues and the Diferences of a Fun^ton, 

15. These relations we propose to investigate on what is 
termed the method of separation of symbols of operation and 
quantity, before applying which it will be necessary to make 
a few remarks by way of explanation. 

By a symbol of operation is meant a symbol put for con* 
ciseness to denote that a defined operation is performed on the 
function, whether general or specified, to which it refers. Thus 
the symbol A used above denotes that two difierent values of 
the variable are substituted in the function and the results 
subtracted from each other. In some cases the function is 
expressed, as where we write Au^; in others we separate the 
symbols A and u^. Thus, in place of 

we write (/>• + PiJ>*'^ +*...+ /?J w^. 

The legitimacy of this step appears from considering the rela* 
tion in which the symbol u^ stands to the symbols within the 
brackets. The symbol u^ is merely a representative of functions 
of a: in general, just as in algebra we denote by a single symbol 
numbers in general : and in the expression Du^ the symbol u^ 
appears in this character, namely as an affix, a frinction for the 
symbol D to operate upon. Hence it may be altogether eiSaced, 
provided we keep in mind that it (or some other representative 



I 

I 
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SEFABATION OF 9YMB01S. 7 

of functions) is that to which the symbol D refers; just as in 
algebra we omit the unit to which a symbol, as (o), refers. 
Therefore also it may be separated from the 2), as in the 

expression 

and this notation is rendered consistent from the analogous 
modes of expression according to which we place the m^ after 
the D in Du^, to denote that the D operates on the w_, and the m 
after the symbol of quantity a in au^ to denote that a multiplies 
the u.. 

By a Symbol of Quantity is meant a symbol put for a 
number or complicated algebraical expression. 

Symbols of quantity may be regarded, however, as symbols 

of operation. Thus if a be a number, it denotes that unity, 

L whether the unit be abstract or concrete, is to be repeated or 

ledded to itself (oj times, and a' or a.a that the same operation 

r IB to be performed on n as was on unity, and o' or o.n' in like 

manner :* and in all cases where the two kinds of symbols 

occur together, as for example in a + A, we must regard a as 

a symbol of operation for the sake of consistency. 

tl6. In the proofs on which are founded a yariety of theorems 
in ordinary algebra, the only properties required to characterise 
the symbols employed are, that they are subject to three laws, 
denominated by a late writer, the commutative, distributive, 
and index laws. 

(1) The commutative, ah = ha, 

(2) The distributive, c{a -^^ V) = ca \ cb. 

(3) The index law, a~.(a') = a'". 



• If in place of the espreasion, miBing of powers, which involves the notion 
of multiplieation, we use as a corroBponding term one implying that the opera- 
tion is performed Buoceasively, it would have this convenience that, we could thua 
the mBBning of the »i in such symbols oa a", i", D", j -j- 1 in the Baine 



8 SEPARATION OF SYMBOLS. 

Hence these proofs hold for all symbols which are subject to 
these laws, whatever be their nature in other respects. Now 
if the symbols which are used in this calculus, A and 2>, (the 
functions they are to operate upon being understood) are so 
subject, then such theorems as those above referred to will hold 
for these symbols. Thus, from the theorem 

(a + bY = a"* + moT'^b + — — ^^ aT-^V 4. . * . 
^ ^ 1.2 

we have also 

(A + DT = A- + mA"-*2) + m.{m •- 1) ^^.^^ ^^ ^ ^ ^ 

1.2 



17. Now that the symbols D and A follow : (1) the com- 
mutative, (2) the distributive, and (3) the index laws, may be 
proved as follows. For the symbol D, 

(1) If ««.,y be a function of x and y, and D refer to re, and 
D' to y only, we have 



D'u = u i« 

(2) If u^y v^ be functions of x, 

(3) As we have before shewn, 

18. For the symbol A, 

(1) If u^^ denote a function of x and y, and A refer to x, 
. wid A' to y,. only. 



GENERAL RELATIONS. 



^^k.y.k - ^..k.y - K.V-* - «..,) = ^X.A., - ^'^'.v = A'A««^^, 

.-. AA' = A'A. 



(2) If «^, t? be functions of x^ 

Aw. + At?. = tt^j + »,,! - (w. + O = A(w. + t?.). 

(3) AX = (2> - 1)X, 

/. A"*AX = (i> - 1 r (i> - !)•«. = (i> - 1)"**X by the proof for D 

= A^^^X, 
/. A*" A" = A*""^". 

19. It might be easily shewn, further, that the symbols — ^ 

2), A, referring to the same function, may be combined in any 
manner conformably with the above laws. Thus, for instance, 

4^ AJD=^DA:=D^A^ &c., 
ax ax ax 

and |(i) + A) = ^2> + ^A. 

20. The definition of D is that Du^ = w^^^, h being generally 
unity. 

Now we know, by Taylor's theorem, that 

d 

where e ** is a symbol of operation referring to w.. Hence 

d 

Also, Aw. = w,^;i - w. = Xfe. - w, - (D - 1) w,, 

. . A = i> - 1. 



10 DIFFERENCES OF ZERO. 

21. To express A"u, by u, and its n successive values. 
Since A = 2> - 1, 

.-. A" - (D - !)• - D" - n.Xr-' + ^l^Lzl) 2?«-' _, &c., 
or introducing a function u^, since 2>"w, = w^^, 

22. If in this series we give to u^ a particular form, we 
easily obtain its n^ difference ; thus let tu = sf*, 

A» (sT) = (;» + n)"" - « (a: + n - 1)« + ^^l (a; + n - 2)"* - &c. 

In this equation, if » > m, A" (a:"*) = 0, whatever be the value 
of X. I{m = n. A" (sf) = 1.2.3. . . .w, and therefore 

1.2.3. ..n « (a; + w> - n(x + n - 1)" + ^'^^"^ ^ (a: + w - 2)*, - &c. 

'^ 1.2 

and if a; = 0, since x may have any value, 

1.2. ..n = «" - w (n - 1)" + ^' ^^ " "^^ (n - 2f - &c. 

V y 1.2 ^ "^ 

Also, making a: = 0, we have 

A^CO"*) « w"* - n.(n - 1)" + ^'^^"" ^^ (w - 2)"* - &c. 

The numbers represented in the form A^.O*", are much used in 
analytical calculations ; they are termed the differences of zero, 
and afford a concise mode of expressing complicated formulae. 

We should find from this formula, making n « 1, 2 . • . . n, 
successively, 

AO"* -= 1, 

AV=2, AV = 6, A'0*=14 

AV = 6, AV=36, AV=150. ... 
where w > m, the value is zero in all cases. 
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23. To express u^^ by w^, and its first n differences. 
Since 2> = 1 + A, 

Z)" = (1 + A)'* *= 1 + wA + ^ A' -I- — ^^ — r-r^ A^.. + A . 

Introducing a function w,, since D^u^ = w,^^ 

«^»+wi - ^« + ^^^« + — ^^ ' ^^* + ..• + wA**'X + ^^«' 

If «^ =i a;"*, u^^ ^ (x ■¥ nTy h being unity, 

.\ (;» + w)^ =t^ (1 + A)" ic", or if a: = 0, 

n"* =^ (I + A)~ O"*. 

a formula by which any number is expressed in terms of the 
differences of zero, 

24. Since A = JD - 1 and D = e***, A being unity, 

A" = (e^ - 1)", 
and introducing the function 

AX = (e^- lyu^, 

which is Lagrange's theorem. 

25. To find the n^ difference of the product of two func- 
tions. 

where D refers to w^, and D to t?^, only. Hence, repeating 
the operation n times, 

1 .iS 









12 herschel's theorem. 

We may expand in a different manner as follows. Since 
Auv^ = (Djy - 1) w.o. = {D (1 + A') - 1} u,v, = (A + DA') uj>^ 

/. AXt?. « (A + 2)A>.t?. 

- {A" + «.A-M}.A' + ?:<!LlLi) A-'l^'.A" + . .} uv 

I 1.2 ^ • ■ 

« t?,A"w, + «At?,A*^*w^j + — —j-;; — ' A'f?,A*^\^, + . . . • 

26. To prove HerschePs theorem ibr expanding any func- 
tion of €^, in a series of ascending powers of t 

Let/(e') = ^ + 5e' + + Cfe^' + 



d 
2 — 



« (-4e^ + Be^. ^ + Ce ''^-^ . .y^,, 



d 



Hence the coefficient of f* is 

fr"* 0"* 

which is HerschePs theorem. 



d 



Cor. Hence the general term in (^ - 1)" u^ is 

1.2. ..m c&* ' 
and therefore 

A '*4y ^ * J. ————— + — — ^— ^^^— ^— — — — — — -4- CZC 

• <fo- 1.2...(n+l) (fo»*' 1.2...(n+2) di"*» 
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which is the general expression for the n*^ difference of a 
function in terms of its differential coefficients. 

27. To express the nth differential coefficient of a function 
in terms of its differences. 

d 

Since ^ = i9 = 1 + A^ 

i = log(l.A), 

or, introducing a function for the symbols to operate upon, 

g=={log(l+A)}-«.. 

This being true for all values of n, writing -n for «, and 
remembering that ( -^ ) is equivalent to I efe". 



/-^«.= {log(l.A)}-X. 



( H ) 



CHAPTER II. 



INVERSE METHOD OF DIFFERENCES. 

28. The Inverse Method of Differences is applied to deter- 
mine from a given difference the function from which it was 
derived. This process is called integration. 

The simplest case is where from the equation 

we have to determine u^ in terms of x. The result is thus 
expressed, 

So that S and A denote operations the reverse of each other. 

Another case is where we have to determine the form of u^ 
from the equation 

AX =/W. 

In this case, the result is written 

u, = S"/(^). 

29. Since any constant quantity connected with the variable 
part of an expression by the signs + or - disappears in taking 
the difference, it must be restored in integrating; and since 
all functions which differ from each other only in the values 
of the constant parts have the same difference, we must, in 
order to include all functions from which the proposed difference 
can have been derived, add to the integral an arbitrary con- 
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stant C, Instead of the constant^ however, we may add any 
function of the variable x^ which remains unchanged when x -^h 
is written for x. This remark will be found of importance. 



30. Since S^A^w, = u, and according to the index law, 
whatever n be, 

A"" A**?/, = AX = '^M y 

2" is identical with A"", and may be interchanged with it. 



31. To find the integral of any rational integral function. 
Since the difference of any rational integral function. 



is a function of the same kind one dimension lower, it appears 
that, conversely, the integral of a function of this description is 
a similar function one dimension higher ; hence, we may assume 
the integral of the above expression to be of the form, 

and then on taking the difference of each side, and equating 
the coefficients of like powers of x, we shall have (w + I) 
simple equations by which to determine the (n + 1) quantities 
PiPi .... /?„+! ; the last term cannot be determined : it is the 
arbitrary constant which is required to give the integral its 
utmost generality. 

Ex. Find S («* + 6x^ + 6a:). 

Let S(a^ + 5ic' + Qx) =p^ 4 p^x + p^ +P^^ ^ P^^y 
therefore taking the difference of each side, 

a:' 4 5a:* + 6a: ^ p^ ^ p^i^^ + 1) +i»8(3a:' 4 3a: 4 1), 

4/?4(4a^4 6a:' 4 4a; 4 1) 
= J»i 4 j»3 4 />3 4 io, 4 (2p, 4 3j»3 4 \p^x 4 (3/?3 4 6^,)a:' 4 4p/, 
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equating coefficients of like powers, 

3/?3 + 6;>, = 5, 3/?3 = 5 - i = |, .-. p^ « g, 

2/>, + 3/>3 + 4;?, = 6, 2;?, = 6 - 1 ^ I = §, /. /?, = |. 

32. Since if w^ = -4^ + -4, a:. 

Therefore taking the integral of both sides, 

whence we deduce a rule for integrating the product of con- 
secutive terms of an Arithmetic Progression, '^ Annex one 
more factor at the beginning, divide by the number of factors 
so increased, and by the common difference.'' 

-n ^/ .r.\/ .\ r ^N (a;- 15)(a;- 10)(a;-5)a;.(a:+5) ^ 

Ex. S(« - 1 0)(a; - 5) xJ^x + 5) = ^ 1-~ ^—^ ^ + C. 

5.5 

It is convenient to resolve by inspection a rational integral 
function into factors as above, in place of integrating by the 
method in Art. 31. 

In integrating by this rule, we must be careful to observe 
that each factor must be derivable from the preceding by 
writing re + 1 for x, thus the expression 

(2a;- 1). 2a:. (2a; + 1), 

could not be integrated by this rule but the expression 
(2a; - l).(2a: + 1). (2a: + 3) could. 
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which furnishes a rule for integrating a fraction, whose denomi- 
nator is the product of consecutive terms of an Arithmetical 
Progression, and numerator constant. It is ^^ Efface the last 
factor, and divide by the number of factors remaining, and 
by the common difference, prefixing a negative sign.** 

34. This method will apply to the case where the numerator 
is a rational function, lower by two dimensions at least than 
the denominator. Thus, suppose the fraction to be 

oar""' + bx"*'^ + .... + Aa: + / 



Ug . , . . «,^., 



assume. 



oaf*"' + hx""'^ + ...4A:ar + /=a' + h'u, + c'uju^^^ + . . . + h'u^u^^^u^^^ . , . u^^_^^ 

Then arranging the second member in terms proceeding 
by powers of x, and equating coefficients of like powers, 
we obtain (» - 1) equations for finding (w - 1) quantities, 
a', Vy c', . . . . k', and the fraction is reduced to the form, 

a' V K 



each term of which is integrable by the preceding rule. 

If the dimension of the numerator were the same as that 
of the denominator, or differed from it only by a unit, we 

should arrive at a term , of which can be integrated only 

approximately. 
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85. To find the integral of a". 

Since Ao- = a**' ^a' = (a- I) a'. 



a' 



.-. So* = + C -. 

a - 1 " ' 

So 2V = ^^ 4 2C , = — ^, + C^.,x + C^., ; 

a - I " ' (a - 1/ " ' " ' 

Similarly S-a» = ^ + ( (7, 4 C;a: + (^.i^^''). 



36. To find the integral of log v^. 
Since 

Alogt?,!?, t?,.j = logt?jt?, t?, - logt?jt?, V, = logt?^. 

.-. 2 log ©. = log t?jt?, t?^_j + log C; 

or as we shall write it, 

S log v^ = log CPt?^.j. 

where Pt?, denotes the product of all the successive values 
of v^, from some fixed term t?,, (or in general t?,) to t?,. 

87. The integrals of trigonometrical functions are readily 
found by substituting their exponential values. 

Thus 2S cos ^e = S {e^'-' 4 .-^^"^ = ^^^ + -^.itzTi ' 

g(«-i)ffv-i _ ^v-1 ^ ^.(.-i)tfv.i _ ^-«#v.i ^ 2cos(a:-l)e--2cosa;g ^ 
" 2 - (e^-' 4 e-^^-^) - 2- 2costl * 

„ ^ cos (a: - 1) 6 - cos :r6 

/. 2 cos xB = \ .^ «r . 

2 (1 - cos 0) 

38. To find the n^ integrals of sin(aar 4 j3), cos (ax 4 /3). 



9 
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Writing - n for w in the results in Art. 10, we have, since 
we have shewn A"** and 2** to be equivalent symbols, 

S- sin (oz + /3) = -?EJ^+4) - J^ OlL^) ; 

(2 sm I ay 



Integration hy Parts, 
39. Since ^^,^, = ^»^^x + ^»m\^^x^ 

which is analogous to the form in the Integral Calculusi 

du 



idxu — ^ uv -\dx 



doi 



For v^ writing St?., the expression becomes 

SwA = «^.St?^ - S (Aw,2t?^^j). 
This formula admits of a very useful extension, as follows : since 

AX^. = {DD' - 1)X^., 
where D refers to w., and D' to t?,, only : we have also, 

AX^ = {^(1 + A') - OX^'. = (A + 2>A')"«A. 

= |a*» + ^^A''■^2>.A' + ^^^^^ A'-^Z)'.A'' + . . . . |«^t?„ 






t?.AX + w.At?.A"-'w,,, + - -S-, AVA'-V., + 



Since this is true for all values of n, suppose n negative. Then, 
since A"" is equivalent to 2", 

SXf. = I'.SX - «Ac.2-"«,„ + ^^^^^ A\S»X« - &c. 

c2 
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Cor. Make n = I, 

2m o = vJ2u - At?^2V,. + AV2V^ - &c. 

t X X M X X*-\ X X-Til 

This formula always enables us to integrate the product 
of two functions, one of which admits of successive integrations, 
and the other leads to zero, as the value of one of its successive 
differences. 



40. To shew that (2> - a)-X = a'-S'^w^o"*. 
Since Da' = a"' = a.a% 

therefore in this case the symbol D is equivalent to the symbol 
a, prefixed to a*. 

Hence, •.* AXa' == {DD - 1)" ua% 

where D refers to a', and D to u^^ only. Hence D' is equivalent 
to a, and we may write, 

AX«' = {Da - \fuX = a'*^ i> - - )\«*. 
T- 1 

Let a = - ; 
a 

then, C^ua'" = a^'^D - afujar". 

Whence, (X> - o)X = a'^^AX^"*- 

Since this is true for general values of w, we may suppose 
it negative, in which case A** is equivalent to 2**, and we have 

if w = 1 , (i> - a)- V^ = a'-'^ua' ; 

a theorem of great importance in the Solution of Linear 
Equations. 

CcUcfuHus of Generating FunMions. 

41. If w, be any function of a:, and the series 
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be denoted by <l>(t); then 0(0 is called the generating function 
of u^f because by its expansion it generates the coefficients 
Wq, Wj .... w^, annexed to the corresponding powers of ^ ; it 
is expressed thus^ 0(Q = Ou^, 



Thus, since log(l -0'* = ^ + - + - + .,,, +-+ .... 

°^ 2 3 X 



.-. log(l -0"'= (?- . 

X 



42. The following theorem is of great use in this branch 
of the subject. 

If 0(Q be the generating function of w^, or f^{t) = Ou^^ 
and i/>(0 any function of t^ which can be expanded in integral 
powers of /, positive or negative, then 

0(0^(0= G^(;^)«^., 
and 0(0 !/'(-)= 04{D)u^. 

Let i/* (0 * .... + «.,^"* + Oq + ^1^ + ^2^ + • • • • + ^.^ + . . . . 
Then, because 

we have for the coefficient of ^, in the product i/'(0.^(0> 
or (..••+ a.jjD + «o + ^1^"^ + flaX)"' + . . . .)«/,, 



or 



(i)"- 



That is to say, Gi/- f -^ j m, = <^ (0-0 (0- = 4' {t)Gu.- 
Similarly, Gi{D)u, = ^p(-\<^(f) = .^f 1") 6?m.. 
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43. A variety of interesting results may hence be deduced. 
Let xP(D) = 2)-", 

then GiD-'^uJ, or Gkc^,^ ^ t\Ou^. 

So if xIj(D)=^I^, 

then 0(D^u,), or Gu^^^ = t'"Gu^. 

Again^ 
Let xP(D)=(D-' I)-. = A". 

then GAX = ( " " m" ^^** 

Writing - n for n, 

Again, let xf, (D) = (log i>)* ; 

then 6?.(log 2))X = (^^87)" ^^' ^ 

but log D = log e^ = -J- , 

Hence GJ^'dx^'u^ = f log - j Gu,. 



Ntmibers ofBemouiUi, 

44. This is an appellation applied to distinguish certain 
coefficients which occur in the expansion of the function 

ajid which are frequently employed in analysis. 
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The development can involve no odd powers of t above 
the first. For, since 



.-. (- =1 ^ + {t)y and .-. (~ - (0 = ^9 

so that vhen in the series for (t) we put - t for t, and subtract, 
all the terms vanish, but that which involves t, which could 
not be if there were any higher odd powers than the first. 

Again, the first term of the series is unity, and the second 
— . For, by Maclaurin's theorem, the first term is (0), and 



2 
since 

t t 1 



^ + _ + ... l-f— +... 



/. 0(0)= 1. 



Again, if (^) = 1 + -4^ + , . . . 

•'• (- = 1 - -^^ + • • • • 
and 0(0-0 (- = 2^^ = -ty :. A = -\, 

Hence the series runs thus, 

_L = 1 _ 1^ + A ^ - J^ ^* + ... + (- \y' ^ f^-\ e- + ... 

tf^-l ^ 1.2 1.2.3.4 ^ 1.2... 2w 

and B^, B^y &c. are called Bernouilli's numbers. 

45. The numerical values of the coefficients may be ob- 
tained as follows. 

By Maclaurin's theorem, if , = u, 

e - \ 






d'w\ t 



3 



2 ^WL 1.2.3 "'•••• 
-j^ ) denoting the value of — , when in it t is put = 0. 
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Hence, therefore. 



^^^^de 



)j "^'=(^lv ^«'($L^ ^"^^^^"' 



and all the odd differential coefficients but the first vanish, 
when in them t is put = 0. 

Now wCe* - 1) = ^ and differentiating n times, 

" 1 when n = 1. 

Putting ^ = 0, and giving n successive values, we have 
w =r 1, w = 1, 

^ , . y 7.6.5.4 , ^ d'u ^ 

»=7, 1-5 + 5- :: • ^ + ^ TiT = <^ > 

^ ^ 1.2.3.4 ^ rf^ 



dt 



» = 9^, -^^ = " iRj> ^^^ so on. 



Hence, ^^ = ^, A = si> ^g = i> ^7 = ij^ 

and any function, which by any means can be brought into 
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the form -j — - , can be expanded in terms involving BernouilU's 
numbers. 



46. To express S^w^ in terms of the differential coefficients, 
and integrals of w^. 

SX = ^"X = (2> - 1)"X = i^^ - i)'X ; 



and therefore by the exponential theorem. 



efe 1.2 dx^ 1.2.3 rfa:' 



j w.- 



47. When » = 1, we can expand in terms of Bernouilli's 
numbers as follows : — 

d 
dx 



s«. = (.- - i)X = ( I"' 



* .«-' 



"•c*^- 1 



,dx) \ Ux'^l.idx' 1.2.3.4 db*"^ '••;"" 



-( 



•^ ^ 1.2 rftr 1.2.3.4 efo' ' 



•^ * ^ ' 1.2 du 1.2.3.4 cfe' 
a series of considerable use in integration. 

48. To prove that if JBj^_i be the n^ of Bernouilli's numbers. 

Since, employing a system of logarithms whose base is e, 

t _ loge' 
7^ " c* - 1 ' 
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therefore by Herschel's theorem, the coefficient of ^ in the 
expansion of this function is 

1 logd + A) „„ 



1.2 2« 



.-/ o". 



Hence, B^., = (- 1 Y" '-^^ O" : 

expanding, we have 

This proposition connects Bernouilli's numbers with the 
differences of zero. 

49. The Differential Calculus is a particular case of the 
Direct Method in Finite Differences, and the Integral Calculus 
a particular case of the Inverse Method in Finite Differences ; 
hence, from results obtained in the latter sciences with an 
indeterminate increment for the principal variable, we may 
pass to results in the former, by making the increment in- 
definitely small. 

50. We shall add a few examples to illustrate this. 

(1) Given A tan a:0 = tan"* j r^^ — j^ , find = . 

^ l'¥(x + h).xff dx 

The problem consists in finding the limit of 

1 . he 



T tan"* 



h l+(a? + A)a?fl'' 

when h is diminished indefijiitely. 

1 . _. M 

1 4 (a; + A) x9 

h0 



Let -J tan* ;; — rr-M = ^> 



.'. tan hu = 



1 +{x + h).x0'' 



tan hu 

,u = 



hu I •¥{x + h) xff 



12* 
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Let h be supposed indefinitely small, 

, ,. . tan hu 

then limit — ; = 1 

hu 

, T . d tan"*a:6l 

and limit u = ;; j7r„ = ; . 

1 + x^ff^ dx 



x^ - hx 



(2) ^x - ^-^ + C, find fi 



It is the limit of SAa;, where h is supposed indefinitely small. 

Hence limit Ythx = limit • + (7 = — + C. 

2 2 

Again, s.:'=_.^^_^(7, 

... 2ic^A=--- A+ ?A^+(7. 

3 2 6 

a:' 
Passing to the limit Jx^dx = - + (7. 

(3) To deduce Taylor's theorem from the formula 

n.(«- 1) , 
^..nk = «^* + ^^«^. + — :[^ ^X + 

This we do by making h evanescent in every term. Let nh = A, 
a finite quantity, then when h becomes evanescent, n is inde- 
finitely increased. 

Also limit of -^ = -=-' , 

h ax 

hmit of -75^ = hmit of ^, ^ = -^-^ , 
A A dixr 

and saon. 
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Now we may write the formula as follows : 

, Aw, nh.(nh - h) AV. 
w,^ - w. + nA -^ + — ^ -^ + 

Hence^ taking the limit of both sides when h is indefinitely 
diminished^ since nh is always. equal to k, we have^ from the 
preceding results^ 

*'* • dz 1.2 da^ 



( 29 ) 



CHAPTER III. 



EQUATIONS OF DIFFERENCES. 



61. An Equation of Differences is a relation between a 
variable, any function of that variable, and the successive values 
or differences of that function. 

Such equations may be thus represented : 

or f{x, u„y Aw„ A'w». . . .A'^Ug) = 0, 

since the forms are convertible; for, as we have seen, u^r is 
expressible in terms of w, and its (r) differences, and vice versd, 
A*'w<r is expressible in terms of u„ and its (r) successive values. 

The order of an equation of differences is determined by 
that of the highest difference or successive value it involves. 

52. The complete integral of an equation of differences of 
the ffi^ order will contain n arbitrary constants. 

This appears from considering the manner in which such an 
equation is formed. 

Let f(x, Ugy a) = 

be an equation between the variable x, the function t^„ and the 

arbitrary constant a. The relation is true when z + h is put 

for X, whence 

f(x + A, w^^^i, a) = 0. 



S2 eouahons of differences. 

It is unnecessary to add a constant after performing the 
integration indicated in the equation 

log Q, = 2 log A,. 
For suppose 

log Q. = 2 log A, + log C = log C'PA,.,. 
Then Q, = CPA,.„ 



u. = CPA.., |: 



and u. = CPA., I S ^r^j- + ^' 



c-l 1 - pA 



^ 



P^ ;s .;=4- + C"(7" 



which is the same result as before, C'C* being equivalent to the 
single constant C. 

Ex. 1. 

Let 

Then, dividing by Q^^j, 



«... ■ 


-e'' 


'tt. = c^. 


e"-'- 




»v 








«. 


~ Q ' 



or A 






But 2a: - 1 = A log Q„ 

/. (a; - 1/ = log Q„ 
the constant being unnecessary. 
Hence Q,„ = c*'. 



and 



(;&)-' 



"' =x+C, 



c(»-'>' 



.-. M, = e*'"' {a: + C). 
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Ex. 2. w^j - au^ = {2x + 1) a*. 



then^ dividing by Q 



'«+i' 



or 







(2a: + 1) a* 
Q 


A 




(2aj + 1) a' 



But A log Q^ = log a, 

• • log Q. = (« - 1) log «f 
no constant being necessary; i.e. 



logQ. 


= log o'-\ 


••• Q... 


= a*. 


A "' 
a-' 


(2x + 1) a* 


o* 



Hence A -^ = ^^ -^^ — = 2a: + 1, 



.-. 3i. = S (2a: + 1) = a:* + c\ 



c' being an arbitrary constant, 

.-. u^ = a*"* (a:* + c'}. 



/fufirdc^ Integrals of Equations of Differences. 

64. We have stated above that equations of differences of 
the first order, when inexact, are formed by the elimination of 
a constant (a) between two equations, as 

K =/(^> ^) ^^^ ^.a =/(^ vhyO) (1). 

D 
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If, however, we suppose {a) not to be a constant, but a function 
of Xy and denote this function by a^, then these equations 
become 

and (since we have made no supposition as to the form of the 
fonction) if we suppose it such that 

fix + A, «^J =/(a? + A, aj (2), 

then the two equations become 

w. =/(ir, aj, u^f^ ^f{x + h,aj (3), 

wherein a^ appears exactly as a did in (1), and therefore, by 
eliminating (aj from equations (3), we shall arrive at the same 
equation of differences as we did by eliminating (a) from (1). 
Now the condition (2) is satisfied by supposing 

«,.;i = ««> or Aa. = 0, 

whence a^ = a, a constant ; and, substituting this value, we have 

but there are other values of a^ which satisfy the condition (2) ; 
for since a^^,^ = a^ is a value of a^^^, which satisfies equations (2), 
^«+* - a, is, by the theory of equations, a factor of that equation, 
and hence we may divide out by it, and the result will be 
(if powers of the constant (a) higher than the first are involved 
in (1),) an equation of differences of the first order with refer- 
ence to a^: by solving which we shall obtain one or more 
values of a^ in terms of x and arbitrary constants, and the 
substitution of these in 

will give other and indirect integrals of the proposed equation 
of differences. 

Indirect integrals of equations of differences are obtained 
on principles very similar to those on which singidar solutions 
of differential equations are obtained, but they differ in this. 
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that whilst the former may contain as many arbitrary constants 
as the complete integral from which they are derived, the latter 
can contain no arbitrary constant. 

Ex. (a: + 1) w^ = xu^, + - (u^,, - uj, 

m 

This may be written thus : 

u = zAu + — (Aw y. 

' m 
Taking the difference, we get 

Au^ «= Aw, + (a: + 1) A\ + - A (Aw J', 

or c= (a: + 1) A'w, + — A (Aw J*. 

This equation is satisfied by 

Aw, = a, 
or w, = aa: -h a'. 

But, from the proposed equation, 

1 2 

w = ax •¥ — a , 

' m 

. 1 2 

.'. a = — a , 
m 

and therefore w^ = aa: + — a^ 

the complete integral with one arbitrary constant. 
For the indirect integrals we must put a, for a. 

Then u = ax + — a/, 

• • m 

.-. w,,, = a,,j (a; + 1) + - (a,J^ 

//• 

d2 
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But by the conditions of such integrals^ 

m 
Subtracting the two preceding equations, 

- (« + 1) (a., - «-) + ^ {(a..)" - («J} ; 

since o^, - a. = is not the solution, 

is the equation for determining the solutions. 

Now a^i + a^ = - m(a: + 1) gives 

- a. = (- If Sm (a? + 1) (- 1)' 

= m?^-(7(-l)-, 
the indirect integral. 



Linear Equations of Differences ofdU Orders. 

55. The general form of such equations is 

^'^ -^ Pi ^^n -^ •••• +/^«««.= X> 

where JO,, /?2.../?^ and X denote functions of a:. To this form we 
may convert the following : 

AX + ?i^"X + + ?„«^. = X* 

where y,, q^* • • •?„ denote functions of x. 

56. We shall first consider the case where the coefficients 
PvP2'*'Pn ^® constants. 
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According to the notation before introduced^ we may write 
the first equation as follows : 

or separating the symbols of operation from those of quantity^ 
we may put it in the form, 

(Z^ + jOjJD"-' + + />J w, = X, 

or, for shortness, f{P) ^* = X. 

Hence, performing on both sides the inverse process oif{I)\ 

«. = f/(^)}-'x = {D- + pjy-' + ....+ pyx. 

Now, let the equation 

^ + />i«""^ + . . . . + A = 0, 

(which is called the auxiliary equation to the proposed), be 
supposed to be made up of a series of unequal factors, such as 
{z - a), and of a series of equal factors such as (2 - J/, where 
(fl) and (J) may be real or imaginary, this comprises every case 
that can occur in the composition of the equation ; then, by 

the theory for the resolution of rational fractions, -^r— may be 

resolved into a series of fractions of which the type is 

M 



9 

z - a 



where 



M I 



dfiz) 



together with a series of fractions of the form 

(z-by'^{z-br'^ "" "^ z-b' 

where N = [-J— -^ ^tlUX 

where iV,, \^,^y^^. ^^^^ ]^^^. 

Introducing the symbol of operation, since the method of 
resolution into partial fractions is independent of any properties 
of the symbol, except that it be subject to the commutative, 
distributive, and index laws, to which laws the symbol D is 
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also subject, we have 

= {Jf(2>-a)-V&c. + [iVr/i)-J)-''+iV^/2>-J)-('->)+&cO + &^^ X 
= {M{D- ayX+&c.}-\-{N,(D-by'X+NlD-b)-^''-'^X+8cc.}'^... 

by the foregoing theorem. 

Each of the integrations will introduce an arbitrary con- 
stant, send the whole function thus introduced, which is called 
the complementary function, will be (making X = 0, and 
replacing the products of constant quantities by single con- 
stants), 

(Aar + ) + (B^ + B^z-¥ + B^)b'' + &c. 

57. The above is the simplest and most symmetrical form 
into which the solution of the equation can be brought, when 
a and b are real. When they are, one or both, imaginary, 
the residt admits of some modification. 

By the theory of equations, if a + /3 V- 1 be an imaginary 
root of the auxiliary equation, there will also be a conjugate 
root a - /3 V- 1. For the sake of generality, suppose this 
root repeated p times in the equation; then the solution will 
consist of a series of pairs of terms, the coefficient of one 
term in each pair being the same function of a + ^ V- 1 
that the other is of a - /3 V-1. By the theory of impossible 
quantities, these coefficients may be represented by (7 + 2> V- 1 , 
and C7- 2>V- 1, so that the corresponding terms introduced 
will be of the form 

(C+ 2>v^-l)(a + j3V- l)-*'S'-(a + j3\/- 1)-X, 

+ (C- 2>v'-l)(a - j3 V-ir'S^(a - /3 V- i)-X, 

and there will be a series of them obtainable by giving r 
all positive integer values up to p. 
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Let a = p COS 0, /3 = p sin 0, 

then (a + /3 V-lT" = p""(co8 + V- 1 sin 0)-', 

= p*-*"{cos(a: - r)0 + V- 1 8in(a: - r)B}, 

consequently the two sums are reduced to 

^^ r[(7cos (a: - r) 61 - Z) sin (x - r) ff] S*" (p-'X cos a:!?) ^ 
^ U [Csin (a; - r) + 2) cos (x - r) B]^"" (p-'X sin xOy ' 

and the complementary function is a series of terms of the 
form 

([Ccob{x - r)0 - 2)sin(a: - r)0](C; + C;:i: + C7^-)l 

^ U [(7sin(a:- r)0 + 2>sin(a:-r)6l](2>o + A^+- + A^r^)^* 

58. The two last Articles exhibit the general solution 
of all linear equations of differences with constant coefficients. 
The process is simpUfied when X in the second member 
has certain forms. 

Thus, if X = 0, the integration is dispensed with, and 
the solution is reduced to the process of finding the com- 
plementary function. 

Again, if X = a^, or the equation is /{D) u, = a**, then 

But since Dd" = a^.a^, the symbol D is equivalent to the 
symbol of quantity a*, which we may therefore interchange 
with 2>, 

a** 
•'• ^x ^ TFT "•" complementary function. 

Again, if X consist of terms of the form Ax"^, that is to 
say, if it be a rational integral function of x^ we must expand 
f{D) or /(I + A) in a series ascending by powers of A, 
so that we shall have 

w, = + . . . , + (^0 + -4iA +.,,.+ -^^A"* + .,..) Ax"^ + • . . . 
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and the series will be terminable, since m being the highest 
index terms after A"*a?~ vanish. 

Ex. I. ^Xf2 ~ 2aw,^j + w. = a^x, 

or (a^iy - 2aD + l)u^ = a^x; 

,'. ( X> — ]u - x; 

:. w^ = [Z> J2: = [A+1 — )a: + comp. fdnction, 

(l — j -2fl — J A+. ...>« + comp. function, 
a^x ?l2a' - 

which function is ( - ] 2'0, or ( - ) (C, + C^x) ; 

Ex. 2. A seed grain is planted, and when one year old 
produces the next year tenfold; and when two years old 
and upwards, produces annually eighteenfold. Every grain 
is planted as soon as it is produced, and the produce of each 
is according to the above law. Find the number of grains 
at the end of the x^ year. 

Let Ug = number at the end of the x^ year. 
Then in the (x + l)^ year there will be produced, 10 each 
by those in their second year, and 1 8 each by all the rest. 



That is, w,,i = w. + 1 (w^.^ - w..,) + 1 8 w, 



«-a' 



or u^,, - t^. - 10 w,.i - 8 w,., = 0. 

Let t^,., = t?,, 

then, t?,,3 - tj^^ - 10 t?^j - 8 tJ. = 0, 

or (iy-Zy*- lOZ)- 8)t?, = 0. 
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The roots of the auxiliary equation are 4,-1, and - 2; 

/. t?.= c;4- + q(-i)'+c,c-2)', 
••• «^. = ^,*a -c,4'+ q(- ly + q(- 2)-, 

changing the constants. 

Making a; => 0, 1, 2 in succession, 

«^i = l = 4(7,- C;- 2(73, 
w^= 11 = 16(7,+ (7, + 40;, 
whence, C, = ^,"(7, = - §, (73 = f, 



and u 



7.4'-27(- l)r + 20(- 2/ 



15 



59. Our investigations hitherto have proceeded on the sup- 
position that the increment is unity : we shall now shew how 
to solve the same class of equations on the hypothesis of an 
indeterminate increment h. 

As before, let the general equation be thus represented, 

then u^ = {f(D)Y'X. 

Let {/{D)}'^ be resolved by the theory for the decom- 
position of rational fractions into a series of terms, 

M(D-a)-'-\- + i\r(i> - 6)-" + , 

then «i, = M(D - ay'X + + N(D - 6)-'X +...., 

a, b, &c. being roots of the auxiliary equation 



42 LINEAR EQUATIONS^ 

Now the theorem from which we find the value of 
JV. (X) - J)"'X is investigated on the supposition of an incre- 
ment unity, and will therefore have to be modified when 
h is taken for the increment. 

We have, as before, 

where D refers to w^, and D' to t?, only. 

Let t?^ c= a* ; then Do, or v^^^ = a*.a*, so that D is equi- 
valent to a*. 

Also AX«?, = {DD' - 1/w.t?, 

becomes, when t?, = a'. 



I 
.A 1 



Let —K- ^9 .'. a = {- ] = a ^ ; 



a* Va 



x+ph 



then (D - a/e^, = a * A" w,o *; 

changing p into - jo, this becomes, 

»-ph g 

(D - a)-X = a^^ S^a"'' . 

Hence, applying this theorem to the general equation, we 
have 

»-h » *-ph » 

w, = + +Ma "^Sa'^X + + Nf^^SPh'^X + , 

the complementary function being of the form 

+ .... + Ifo* +....+iV^J* f Co+CjT + . . .. j. 
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or + JfV+ .... +N'b^(c,-^ C,^+ ...\ 

changing the constants. 

This may be written 

(ly - 2mD + m* + n*)t<, = ; 

Cj and (7j being arbitrary constants. 

Let m = /o cos 0, n = p sin ; 

* MX X 

:. u^ « Cj/o* (cos + V- 1 sin 0/ + Cy (cos - V- 1 sin 0/. 

= /t>^{(C.+ C7,)cos ?0 + V-l (C, - C7Jsin|0J, 

= p* I -4 cos -7- + -o sm — J , 
changing the constants. 

60. The two following examples are added with a view 
to illustrate the analogy that exists between the methods and 
results in the Differential Calculus, and those in the Calculus 
of Finite Differences. 

Ex. From the fundamental theorem 

X X 

to pass to the corresponding result in the Differential Calculus. 
We have, since 2> = 1 + A, 

x\ X 

{A - (a - \)Y^u, = a'^'^YFu/^, 
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Let a - 1 = oA, 

.-. (A - oA)-' = (1 + oA/"' 2^W,(1 + ah)"\ 



■■• [t - -r 



(l+aA/"'s''AX(l + «A)"*. 



Passing to the limiting case^ when h is indefinitely small^ 

since Umit of -r- = -- ; and limit of S'^.A'' = j daf ; 

h ax J 



~h >»-«* 



and limit of (1 + ahf = (f; and limit of (1 + oA)"* = e 



61. A Differential Equation is a particular case of an equa- 
tion of Finite Differences, viz. the case where the inde- 
terminate increment is indefinitely small ; and from the solution 
of the latter, we may pass to the solution of the corresponding 
differential equation. The equation whose solution has just 
been obtained wiU furnish an instance. 



This equation is, when expressed in terms of At^,, A\ 
AX - 2 (w - 1) Aw, + {(m - 1/ 4 »'} w, = ; 

Let w - 1 = m'hy and n = r^h, 
then ^ - 2m' ^ + (m'V n") u^ = 0. 



sy 



Passing to the limit, when A is indefinitely diminished, 

dx dx 
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The solution is^ 



u, = p^lA COS -rr- + -B sin — J : 



but p' = m» + »' = (I + m'hf + (n'hy ; 



,'. p^ = {(1 + m'hf + (»'A/}^ 

tberefwre limit of p* = limit of {1 + 2m' h + (m'* + n'*) A'}* 

-limit of {1 + 2m'A}^ 



c"**; 



and cos -7- = cos 7 tan * — = cos hr tan 7= ) , 

therefore limit of cos -=- = cos tIx\ so limit of sin -^ = sin rix ; 

h h 

.'. u, = e"^' {A cos n'x + 5 sin n'x)» 



62. No general method has hitherto appeared for solving 
linear equations of diflferences with variable coefficients. In 
certain cases they can be reduced to equations with constant 
coefficients^ in others they can be solved by finding the integral 
of an equation derived from the proposed, by destroying or 
modifying certain of its terms, and completing the solution 
by variation of parameters, as in diflferential equations, and 
a large class of equations be solved by the aid of the following 
theorem. 

68. If a particular integral of a linear equation of any order, 
with no term independent of u^, be known, the equation may 
be reduced to one of the same kind of the order immediately 
inferior. 



46 REDUCTION OF EQUATIONS. 

Let the equation be reduced to the form 

A-w. + p^^''-\ + + pjii^ = 0, 

and let u^ = t?^ be a particular integral. 

Assume t/. = t? ^w . 

Then A** (t?.SM?J = (DD' - 1)» (t?,SM?J, 

where Z> refers to t?„ and 27 to St£?., only. 
Hence AXv^^wJ = {2> (1 + A') - 1 }X2tt?, = (A -f 2>A')**t?,Sw?. 

= {A" + nA*-^2>A' + "^-^^"-^ A"-^2>'A'» +....+ 2>»A'"} « Sw; 
^ 1.2 J • . 

= A"tJ,St(?^ + n/y^^^Ws + ^'^^^ A""'t?,^2At(?, + ... + t?,^A""'tt7,. 

1 .^ 

For » writing (n - 1), (« - 2)...&c., in succession, and sub- 
stituting the corresponding values of A**u„ A**"*w,, &c., in the 
proposed, it becomes 

Tl (ft "— O 

A"t?,SM?, + wA*'"*t?,^jW?,+ — ^ A"'*t?,^At(?, + &c. 

+ j!?i{A*''^t?,2w?,+ (n - 1) A""*t?,^jW?, + &c. 

+ (w - 1) At?,^.,A'"'t(?, + t?,,^.^A*-'t(?,} + &c. 

Or (AX + PA'*^ + &c. 4 p^v,) St£?, 

f n(n-l) , (n-l)(n-2) a. ^ ., \ a./, 

+ I ^ 2 '*' r2 '"» + A-a^'W ^^^'^ 

+ &C. 



VARIATION OF PARAMETERS. 47 

Or^ since the first line is zero^ we get an equation of the form 

A**'^w, + yjA^^'tf?, + &c. + y^.jW?. = 0, 

?i' ?3' ^^* heing functions of x ; and this equation is of an order 
immediately inferior to that of the proposed^ and of the same 
kind. 



64. CoR. Since w. « t?.2t(?., .*. «?„ = A — : 

V 

m 

and if we know another particular integral of the proposed^ 

u = vj, then 

/ 

t(?^ = A -^ 

is a particular integral of the last equation, which may be 
reduced now one order lower, that is to say, the solution of the 
original equation may be reduced to that of one, two orders 
lower: and, generally, if we know r particular integrals of 
a linear equation of the n^ order, its solution may be reduced 
to one of the (n - r)^ order. 



Variation of Parameters. 

66. This method, which is extensively employed in the 
solution of differential equations, may be applied also to equa- 
tions of differences, as in the following example. 



To solve w,., - A u^ = J5.. 



jr+l 



The solution of u,^^ - Aju, = 0, 

or "^ = A,, 

is w. = CPA,,,, 

where C? is a parameter. 
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Now suppose Cj in place of being constant^ is a function 
of Xy denoted by (7,^ and assume that the integral of 

u , — A.U ^ S 

is of the same form as that of the modified equation ; that is^ 
assume the integral to be 

where the form of C^ is to be determined from the condition 
of its being the solution. 

Then «,,, = C,,,PA., 

and «„, - A,u, = PA, (C,„ - CJ, 



' PA ' 

m 

and .*. (7 = S ( -zpr-r- + C' | , C a constant. 

Whence, u, = PA^, |s -^ + c[ . 



SimtUtaneom Equations, 

66. In these, instead of one equation of differences from 
which to determine a single function of Xy we have two or more 
equations, each involving two or more ftmctions of Xy and their 
differences or successive values. 

Where the equations of this class are linear with constant 
coefficients, the principles developed in the preceding pages can 
be applied to their solution. The theory of the method is as 
follows : since the symbols of operation employed in this Cal- 
culus are subject to the same laws of combination as those of 
quantity, they may be treated in the same manner. Hence we 
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have only to separate the symbol of operation from its subject^ 
and then proceed to eliminate as many as possible of the 
functions from between the given equations^ treating the 
symbols of operation exactly as if they were symbols of 
quantity. The difficulty of elimination thus becomes reduced 
to that between ordinary and algebraical equations. 

The following are illustrations of the method. 
Ex. 1. w,^j = 2t?. + 12w^, 

These may be thus written : 

(ly - 62> - 14) v^ - iDu, = 0. 

Eliminating w^, regarding the coefficients as if they were 
symbols of quantity, 

{(2y-62>- 14) (2>- 12)- 142>}t?. «0; 

reducing and simplifying, 

(Zy - ISjy + 442> + 168) f?^ = 0, 

or {D - 6)(2> + 2) (2> - 14) t?. = 0. 

Whence t?. = a6* + i (- 2)' + cCU)*, 

and u, = — t?. 

* - 2a6' + f J(- 2y + c(14)**^ 
Ex. 2. w.,, 4 2t?.^, - 8w, = a«, 

£ 
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Separating the symbols, 

(ly - 8) w, + 2Dv^ = a', 

2)e^ (2>» ^ 2) « = - i . 
Eliminating t?,, 

{(jy - 8) (Zy* - 2) + 22>»} w, = (2y-2) a'- 22> i , 
thatis, (i)'-27 w. = (i)'-2)a'- 21)1 . 



Whence, since Da' = a.a*. 



2 1 



w = y-^ — —. <f - ,^ .A + comp. function. 

Now 
(ly - 2'r = ife(-D+ 2)-' + i^(i> - 2)-' + ;fe(J5 + 2)-' - :^(i> - 2)"', 
•. {If- 2ro =^(C;+ C7,a:)(- 2)- +^(C,+ C^)2'+:^C,(-2)'-^C.2-, 

is the complementary function. 



Also « = 



a--(i)^-8)t^, 
2D 



The number of constants, though apparently six, it will 
be seen on inspection, can be combined so as to give only 
four. 



Application of the Calculus of Generating Functions to the Solution 

of Equations in Finite Differences. 

67. The Calculus of Generating Functions has been applied 
to the solution of linear equations of diflferences. The method 
of solution is founded on the following theorem, demonstrated 
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in a former article : 



G^(iD)u, = rp(l\Gu,. 



Now let/{D) t<^ =3 X be an equation of differences : then 

Let F(t) be the generating function of X ; then, applying 
the previous theorem to the first member of the equation. 

Let C^ be the coefficient of t in the expansion of . Then 

Urn" C,\a the solution of the proposed. 
Ex. Let the equation be 

or (D* - 3iy + 3X> - 1) M. = a; j 

i.e. {D-\fu, = x, 

:. G(D-lJa,= Gx, 



or 



l-\\Gu, = t.{\-ty'i 



■■■ Gu. = i^(i-ty' 






£ 2 



64 PARTIAL DIFFERENCES. 

On Partial Differences. 

70. If u denote a function of two variables x and v. the 
expression 

is termed the total difference of the function u : and the 
expressions 

are termed the partial differences of u^^ with respect to x and y 
respectively. 

In like manner, the expression u^^^j^ is termed the first 
total successive value of w,^, and the expressions «<«*A.y»^»y+* ^^® 
termed the first partial successive values of u^^ with reference 
to X and y respectively. 

We shall use the following notation : 

Also 2^, 2y for summations with respect to x only, and y only, 
respectively. 

The values of h and A are generally unity, unless the con- 
trary be expressed. Since 

d d 

therefore D is equivalent to e ^' ^^ , 

d d 

h — k — 

Also D^, Dy are equivalent to e ^' and e *'*' respectively, 
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71. The symbols here introduced are subject to the three 
fundamental laws of algebraic combination ; and hence a variety 
of relations between partial differences and successive values 
might be deduced, as in the case of ordinary differences. 

72, It is obvious also that, as 

when D and A are total operations referring to both variables, 

and, similarly, 

y. ^ fi,(n - \) .^ 



Equations of Partial Differences. 

73. An Equation of Partial Differences is a relation between 
a function, its partial successive values or differences, and the 
variables or constants involved in the function. 

We shall consider only Linear Equations with constant 
coefficients, and such as may be reduced to these. 

By the method of Separation of Symbols the integration of 
linear equations of partial differences with constant coefficients 
is reduced to the same processes as apply to ordinary equations 
of differences of the same class. 

There is an important distinction, however, to be observed 
with respect to the parameters in the two cases. As in the 
solution of ordinary equations of differences we continually 
meet with complementary functions of the form 
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80 in solving Equations of Partial Differences we shall find 
expressions of the form 

in which the arbitrary function (7, takes the place of the 
arbitrary constant C. 

Now we have before shewn that 

Hence, in the solution of Equations of Partial Differences^ 
arbitrary functions of binomials play the same parts as arbitrary 
constants multiplied by exponentials do in ordinary equations 
of differences. 

Ex. 1. W..1., -«/.„,! = :i: + y. 
This may be written, 

(D, - D^) u^^^ = iP + y, 

whence, w..^ = B;^^JD^{x + y). 

The integration implied in the S« will bring out not a 
constant, but a function of y, for the summation in S^ extends 
only to X, 

Hence u^ ^ = ^/"'S,y = D/"* {xy + (7 ), 

= a;(y + a:- 1)+ C;,,.^, 
«a;(y + a;- 1)+ CT,,., 
where we change (7^^.j into C^^ for the sake of elegance. 

Ex.2. t/..2„ - 2w,,j,^,, + t/,,^,, = a;. 
This may be written, 

or, (A - A)X. . = «• 
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Whence u ^ = D/^ S '.2> "x, 



Jx-2)(x-l)x 



changing the functions. 

Ex.3. u^^^-a\^ = 0. 
This may be written, 

(D,'-a'D')u =0. 
Whence, 

«-„ = J A" {(D, 4 aD,)- + (D, - aDX} (0), 
= iD.-' {(- aD^rC, + (aD,)-C,}, 

- U(- «r <?.«-,+ a-'C^«}, 

changing the functions. 

This may be written, 

or (D,-a)(D,-6)«„, = C-'. 

Whence, «,., - (D. - a)-' (X>^ - J)rV», 

= ^ _^ ^ + comp. function ; 

(c - a){c - 0) 
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(since Djf"^ = c.ff"^ = Djf^, and therefore D^ and D^ are each 
equivalent to c), 

•" (c-a)(c-i) " 

Ex. 6. «. , - w s= a*"*. 
This may be written, 

or ^;(i>.-l>r)t',., = «-^; 






^r yfn ^l-« _ | "*" ^y-fi«+i»» 






This result might have been obtained at once by remem- 
bering that 2>, is equivalent to a, and D^ to o"\ 



74. The foregoing principles will enable us to effect the 
solution of the more general equation, 

where the coefficients are constants, and V denotes a function 
of X and y. 

The equation may be written as under, 
(A" ^pPr'D, + P^DrD,' +....+ pD«)u,, = F. 
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Whence, 

The auxiliary equation, 

^ ^ p^'^ + +/>„= 0* 

will in general be composed of unequal factors, real or imaginary, 
one of which may be represented by z - a, and of equal 
factors, real or imaginary, one series of which may be repre- 
sented by {z - by. 

Therefore, introducing the symbols of operation, we shall 
have 

M..^ = {( Ai>;' -«)... .(2).i)/' -bf.... }-'2>,- V, 

and by the theory of partial fractions this may be put in 
the form, 



M,„ = {M(D,D,' - a)- + . . . . 
= M{D, - ai),)-'2)/'"'' r+ .... 
Hence, 

The complementary functions are supposed to be included 
under the signs of summation; to see their forms we need 
only write zero after those signs. 
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Therefore complementary iunction, 
- M(aD,y-'2X0) + 

+ iV;(J2),r''S.'(0) + .... 

+ NXhDj'(,C^ + ^c;a: + . . . . + ,..C/-), 

T • • • • 
I » • • • 

or changing the form of the functions, 

+ N^lf-'QC^,, + ,(7,,^ + . . . . + ,.,C^„a;'-'), 

I • • • • 

If the quantities corresponding to a ... i ... be imaginary^ 
the results will have to be modified, as we have shewn 
in the case of ordinary equations of differences. 



Of the Failure of FormtdcB in the Inverse Processes of the Calcfdus 

of Finite Differences. 

75. There are several instances in the Calculus of Finite 
Differences, as in the Differential Calculus in which a certain 
relation of the subject to an inverse operation makes the 
general formulae fail. The method of arriving at the true 
value of the integral in these cases of failure is as follows: 



CASES OF FAILURE. 61 

since the function in this particular case becomes infinite^ 
we so assume the arbitrary constant in the complementary 
function^ as to make the formula for that particular value 
which gives rise to the discrepancy take the form g, the true 
value of which must be found by the usual rules. The con- 
stant is thus made negative and infinite^ so that the difference 
of two infinite quantities may be finite. 

We proceed to consider the two following examples : 



(I* 
Ex. 1. The expression 2a* «= + c fails when a = 1 ; 

required the true value of the integral. 
Assume c' - c - 



a- 1 



then So* = + c, 

a - I 

Now limit of when a = 1 is same as limit of -^- — , 

a - 1 1 

or its value is x ; 

:. Sa* = a: + c when a = 1. 



Ex. 2. The solution of the equation^ 



18 w. = T -i-car, 

o -^ a 



M 



which fails when J » a required the true value of u^. 

Assume, c' » c -- -r , 

— a 
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then u = + CO*. 

' J - a 



= X(f'^ + CO*, 

when i = a. 



Ex. 3. tt^^j - 2aw^j + a'w, = J*. 



The solution is 



J* 



which fails for i =» a required the true value of u^. 
T '- 1 

umit of — r„ = limit of — -; r = limit of — ^^ , 

{b-af 2{b-a) 2 ' 

2 ' 



Functwncil Equations, 

76. There is a large class of functional equations of which 
particular^ and in some cases general^ solutions can be obtained 
by the aid of the Calculus of Finite Differences. The following 
examples are added with a view to illustrate the method 
of applying it for this purpose. 
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Ex. 1 . F(a + x).F(a - a;) = a' - a:*. 

To find a particular solution. 
This may be written, 

F(a + x) F{a - x) ^ ^ 
a -\- X ' a " X 

or f(a + x)/(a - a:) = 1, 

•r F(a + x) r,. . 

a + x '^^ ^ 

Let a + a: = w^^j, 

a - a; = w^, 

.-. (jt> + l)w, = 2a. 
Hence w, = a - (7,(- 1)" = a - a:, 

where C7 is a function of z, such that 

(7, may therefore be an arbitrary function of sin 2Trz and 
cos 2Trz; to get, however, a particular solution, suppose C^ 
a constant, and replace it by the letter C, 

:. X = C(r 1)*. 
Again, let /(a + x) = t?,,„ 

f(a - a;) = t?., 

••• log «^.fi + log ^. = 0* 
/. (2>+l)logt?.= 0, 
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••• log t). = C" (-!)•= ^ X, 

e mxy suppose 

/. fix) = ^^'"'\ 
a result which may easily be verified. 



Ex.2. / (i±£j - a/(ar) = 0. 



Let X ^ u and = w.^,. 

■ \ - X 



Then eliminating o;^ and multiplying up. 



^..v^. - «*..! + w + 1 = 0. 



The solution of this is 

a; or tt, « tan (C + \wz), 
C being an arbitrary constant. 

4 
Whence ;? = - (tan"^ x - (7). 

TT 

Again let f(x) = t?,, 
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* -1 

or fix) = cue 



tan ^a 



changing the constant. 

Ex. 3. Beqidred a particular solution of the equation, 

(a:) + (a; + (x)} = 0. 
Let «{^+ «(^)} = w.,i, 

and C^) = w, 5 



the solution is. 



u, = c(- l)' = «(a:), 



9 



where c may be a function of 2;, which remains unaltered 
when ;? + 1 is put for z. 

Again, let x + i> (x) = t?,^„ 

X = V 

.-. 0(a:) = o(-l)' = t?^,j-t?,, 
whence, «?, = — + a = a?, 

£; and a being supposed constants, 

/. i>(x) = c(-iy = 2(a-x), 

a particular solution. 



( 66 ) 



CHAPTER IV. 



ON SERIES. 

Inte^ation of the General Term, 

77. A series is defined to be a regular progression of terms, 
increasing or decreasing in magnitude according to a certain 
law. 

The position of any term in a series is distinguished by the 
number of terms by which it is removed from some particular 
term selected for the purpose. This number is called the index 
of that term. 

Hence when the law of formation of a series is given, and 
the index of any term assigned, the magnitude of that term can 
be computed, and thus the successive terms of the series pro- 
duced in order. 

78. Where we wish to regard a term with reference only to 
the manner in which it is produced, and not with reference to its 
numerical value, we put for the index a general symbol, as x, 
and the term thus exhibited is called the a^y or general term 
of the series. 

79. The law of formation of a series is represented by the 
function, which expresses the general term in terms of its 
index. Since the magnitude of every term depends only on 
its index, and the law of formation of the series, it follows that 
the form of this function does not change in passing from one 
term to another. 
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80. Wherever the general term of a series is an integrable 
function^ we are enabled to find the sum of any number of 
terms taken in order, by means of the Calculus of Finite 
Differences. 

There are two cases to consider: first, where the general 
term is given as an explicit function of the index ; and secondly, 
where it is given in terms of the index by an equation of 
differences. 

81. Case 1. Let /Si be the sum of x terms of a series, whose 
general term is u^. 

Then 5 = w, + w, + ....+ w , 

«, denoting the first term. 

Also, -S^,, = w, + Wj + + u, + M„,. 

Hence S,^, -S, = m,., = AS„ 

■■■ S^ = Sm,„ + C; 

C denoting an arbitrary constant. Hence the sum of any 
number of terms ending with the general term w^, is equal 
to the integral of the succeeding term together with a constant. 

To determine the constant, make a; = 0, 

or = SWj 4- <7, as we shall write it, 

••• ^. = Sw,,, - Sw, ; 

and in general the arbitrary constant will be determined by 
the term with which the series commences. 

Ex. 1. To find the sum of x terms of the series 

2+12 + 36 + 80 +&C. 

f2 
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The general term is 

• 4 3 

TT « (x-l)x(x-^l)(x-¥2) x(x-^ \)(x-\-2) ^ 
Hence 2w^, = ^^ ^ ^ , — — + — ^^ ^ + C, 

4 o 

.-. 2w, = (7. 

Hence S = (^"O^C^-^ l)(^+2) , ^(^+ l)(a: + 2) 
■ 4 3 

Ex. 2. To find the sum of x terms of the series 

1.3.5 + 3.5.7 + 5.7.9 + . . .. 
Here u^ = (2x - 1) (2x + 1) (2x + 3), 

" ■* 4.2 ' 

_ (2x - 1) (2x + 1) (2a: + 3) (2x + 5) 
•*• ^^»*i = > 

... sum to:, terms =.^^^^^)^^^-^^)^^^^^^)^^^^^^)^^^ 

8 

Ex. 3. To find the sum of x terms of the series 

1.3 3.5 5.7 
Here w, = 



• (2a;-l)(2a: + l)' 
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1 

2(2a;-l)' 

Sw-*, = ; ;^ ai^d Sw, = - i 

•*' 2(2a? + l) ' 

Therefore sum to x terms « A z : , 

* 2(2a; + l) 

and sum to oo » ^. 

Exk 4. To sum a; terms of the series 



1 3 5 

+ + + . . 



1.2.3.4 2.3.4.5 3.4.5.6 



Here u. 



' a; (a: + 1) (a; + 2) {x -v 3) 



1 



{x + 1) (a; + 2){x + 3) a; (a; + 1) (a; + 2) (a? + 3)' 



_ 1 2 1 1 ^ 

•*• • " 2 * (a; + 1) (a; + 2) "^ 3 • a:(a; + 1) (a: + 2) "^ ' 



1 1 1 ^ 

*'• •*>""3*(a;+l)(a: + 2)(a: + 3) (a: 4 2) (a: + 3) "^ ^ 



and 2w, = - . ----- - r-r + ^« 

* 3 1.2.3 2.3 



Therefore sum of x terms 



3\(a;+ l)(a; + 2)(a; + 3) 1.2.3J l(a: + 2) (a: + 3) 2.3j' 



2 1 1 
and sum to oo = -- . — = - 

3 2.3 9 
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Ex. 5. To find the sum of x terms of the series 
a cos + 2a' cos 20 + 3a^ cos 30 + ... . 
Here u^ = xa' cos a:0, 

and .-. 2w^ = arSa* cos xO - 2V cos (a; + 1) 

a cos (a: - 1) - cos a:0\ 



= a: I a*. 



a^ - 2a cos 



~^^ ) 



_ a*, a*'* cos (a: - 1) - 2o*. a' cos a:© + a*"^^ cos (a: + l) ^ ^ 

(a* - 2a cos + If 

^ / , \ «fi ^ cos xd - cos (a; 4 1) 

.-. 2w,^. = (a: + 1) a''\ , \ ^ 

'*' ^ ^ a' - 2a cos + 1 

^^2 a^ cos a:0 - 2a cos (a: + 1) + cos (a: + 2) -^ 
* (a*- 2a cos 0+ 1/ **" ' 

_ a - cos B 2 a' - 2a cos + cos 20 ^ 

^ a' - 2a cos + 1 (a* - 2a cos + 1/ 

(iJ? + 1) a*^* {« cos a:0 - cos (ar + 1) 0} - a (a - cos 0) 
'" a* - 2a cos 0+1 

a* '*{a'cos a;0-2a cos(a:+l)0+cos(a:+2) 0}-a^a*- 2a cos 0+cos 2fl) 

(a'- 2a cos + lf ~ ' 



Recurring Series, 

82. Case ii. In this case the first step towards the sum^^ 
mation of the series is to solve the equation of diflTerences 
which gives the form of the general term: to do this in 
all cases is not possible ; where, however the equation is linear 
with constant coefficients, the solution can, as we have seen, 
be effected. Series corresponding to this last class of equatians 
jwre called Recurring Series. 
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83. The equation of differences which determines the 
general term is in this case comprehended in the form, 

and the series of coefficients which connects any term with 
those which precede it, is called the Scale of Relation. 

Thus 1 +ji?i + .... 4^?^, 

is the scale of relation of the series, the general term of which 
is given by the equation written above. 

84. If a series be stated to be recurring, the number n 
of terms in its scale of relation being given, the scale itself 
must be found by assuming an equation of connection, 

and by giving to w,, w^^j, .... the values corresponding to the 
given terms, so as to form a member of equations sufficient to 
determine the coefficients p^^p^ ' ^ ' Pn' 

85. A recurring series may in most cases be considered 
as made up of a series of Geometrical Progressions. 

For if a,, ^2 ... . a^, be the roots of the equation, 

2" + p^s^-^ 4 .... +p^= Oy 
which is auxiliary to the equation of connection, 

and these roots be all real and different, the solution of the 
equation last written is 

where Cj, c, . . . . c^, are arbitrary constants. 

Hence, making a; = 1, 2, 3 ... . successively, and collecting 
the terms multiplied by the same constants, the series 

^1 + ^2 + • • • • + W, + .... 
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is transformed into 

<?, («! + a,' + (jj* + .... + a'\ 
+ c, (o, 4- o/ + 0/ 4- + o/), 

"1 • . • • 

which is made up of a series of geometrical progressions. 

86. Where however the auxiliary equation contains equal 
or impossible roots, the form of solution is altered, and the 
proposition no longer holds. 

For suppose the equation to have r roots equal to Oj. 
Then w, = {c^ + c^x + cj^ + ... + c^.iOf" *) a^ + c^+i«'r+i + •••> 

and on making a; = 1, 2, 3 ... . successively, we do not resolve 
the series into progressions as above. 

If there be imaginary roots, the solution takes the form, 

«, = C^o + CjO; +. . . • + ^r-i^*) P* cos x9 

+ (c/ + c/a: + . . . . + c'^.ja;*"^)p*sina;0 +....+ c^a/. 

In this case therefore the series cannot be resolved into 
geometrical progressions. 

87. The solution of the connecting equation will involve 
a number of arbitrary constants, to determine which, the general 
term must be made to coincide successively with a number 
of the given terms. 

Ex. 1. To find the sum of x terms of the series, 

2 + 16 + 68 + 244 + 

^hose scale of relation is 1 + /?j + /?,. 
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Let the equation be, 

Hence, 68 + 16/?^+ 2/?, «= 0, 

244 + 6Sp^ + \6p.^ = 0, 
from the fitst^ 

544 + 128/?j+ 16/?, = 0; 

.-. 300 + 60/?, = 0, .-. /?! = - 5, 

and 2p, - 80 + 68 = 0, /. p, = 6, 

and «,^, - 5w,^, + 6t^, = ; 

or (D* - 52> + 6) t^, = 0, 

or (2>- 3)(2>- 2)t/,= 0, 

w, «= c,3* + c,2*. 

Also w, = 2 = 3c, + 2 Cj, 

^2= 16 = 9c, + 40,, 
whence ^i = 4, c, = - 5. 

Hence, t^, = 4.3* - 5.2*, 

3* 2* 
/. Sw, = 4. 5. — , 

2 1 ' 



2w,,, = 2.3*"' - 5.2**', 

2t^, = 2.3 - 5.2 ; 
.'. 1?,= 6.(3*- 1)- 10(2*- 1), 



the sum of x terms. 
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Ex. 2. To sum z terms of the series, 

1 + 6 + 29 + 132 + . . .. 
whose scale of relation is 1 +/>, +/?,. 

Let w,,, + PjW,,j + p^u, = 0, 

be the equation. 

Then 29 + 6pj + p^ = 0, 

132 + 29p, + 6p, = 0. 

Multiplying the first equation by 6, 

1 74 4 36/?, + 6/?, = 0, 

whence, 42 + 7pj = 0, /. p, = - 6, 

and Pj, - 36 + 29 = 0, /. p, = 7 ; 

or (I?- 62>+ 7)w^=0. 

The solution of jD' - 62> + 7 = 0, 

or (Z) - 3)' - 9 - 7 = 2, 

is 2> = 3 + V2, and 2> = 3 - V2. 

Hence, u^ - c, (3 + V2)' -h c^(3 - \/2)*. 

Making a: = 1, 2, successively, 

«, = 1 = 3(c, + c,) + V2 (c, - cj, 

w, = 6 = 11 (c, + c,) + 6V2(c, - c,), 
from the first. 



.*. = Cj + c^. 
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1 1 

=^ C, - C„ .-. Cj =• - C, = 



/. 2V2 W, = (3 + V2/ - (3 - y/2y, 

.-. 2V2 S.„. = (^- ^^2^ - (« - ^T > 
'^' 2 + V2 2 - V2 ' 

„ /- ^ 3 + V2 3 - V2 

2v2 Sw, = 7- , 

' 2 + V2 2 - V2 ' 

*"* » 2V2 \ 2 + \/2 2 - V2 J 



Convergency and Divergency of Series, 

88. A series is said to be convergent, when any number 
of its terms being taken, the sum of these terms has a finite 
quantity for its limit, when the number of terms so taken 
is indefinitely increased ; and divergent in the contrary case. 

89. In a convergent series, the sum of any number n 
of consecutive terms tends continually to zero as n is increased 
without limit. 

Suppose 5^ = Wj + ^^2 + .... + ^^^. 

Then, «n^ = «*, + «^2 + ••••+ «^„ + K.i +....+ u^^ 

and /. 8^^^ - «, = u^^^ ^u^^^-^- .... -^^ u^^^. 

Now, as w is increased indefinitely, each of the quantities 
8^^^ and s^ approach to a finite limit, which limit is the 
same for both. Therefore s^^ - $^ approaches to zero as 
its limit ; which proves the proposition. 

Making m = 1, we see that u^^^ continually tends to zero 
as n is increased, which however is not a sufficient test of 
convergency. 
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Tests of Convergenot/, 
90. (I) If w be the n^ term of the series, 

Wj + Wj + .... + w^ + . . . . , 

and the ratio -^^ approach to a finite quantity / as its limit, 

when n is indefinitely increased, the series will be convergent 
or divergent, according as / is less or greater than unity. 

Let i be a number less than the difierence between 1 and /, 
so that the two quantities I - i and I + i may be both less 
or both greater than 1 , according as / is. 

Then, by assigning to w a value large enough, the ratio 
-^ may be made to lie between I - i and / + t for all superior 

n 

values of n, and the terms of the series. 



«^«» K.x> K.2y &c., 
will be confined within the corresponding terms of the series, 

and u^, u^(l + t), u^ (/ + if, &c. 

Now both the series last written, and therefore the in- 
termediate series, 

K> K.V «^«.2» &c-> 

will decrease indefinitely if u^ does so, as n is increased without 
limit, supposing / is less than 1, therefore the series is con- 
vergent in that case. 

Also both series, and therefore the intermediate series will 
increase indefinitely if I is greater than 1, or the series is 
divergent. 

If Z = 1 , the test gives no result. 
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9L (2) The series, 

is convergent, or will become so, if the superior limit of 
1 

(wj* be less than unity, when n is increased indefinitely, and 

divergent in the contrary case. 

Let L denote the superior limit of (w^y* when n is infinite. 

Suppose first L less than 1, and let A be a number between 

1 

L and 1. Then, on increasing n without limit, (uj* will 
approach so near to L, as to become constantly less than A, 

Hence we can give n so large a value, that for it and 

1 

all superior values, (uj** is less than A, and therefore u^ less 
than A" ; consequently the series will conclude by always having 
its terms less than the corresponding terms of the progression, 

which is convergent, since k is less than 1 ; therefore, d fortiori^ 
the series proposed will end by becoming convergent. 

In the next place, suppose L greater than 1. Let k be 
a quantity between 1 and i, so that h is greater than A, 
and k greater than 1. Then, when n is increased indefinitely, 

(t^Jf on approaching indefinitely near to h will become con- 
stantly greater than A. Hence we can give n so large a value, 

that for it and all superior values we shall have {uj" greater 
than ky and therefore u^ greater than A**, and the proposed 
series will conclude by having its terms greater than the 
corresponding terms of the geometric progression, 

which is divergent, k being greater than unity. Hence the 
proposed series will end by being divergent. 

\ili- 1, this test also gives no result. 



78 CONVERGENCY AND DIVERGENCY. 

92. The limits / and L determined in the two last Articles 
are identical^ i.e. 

limit oi(^^ = limit of (w J».. . 
For, hy taking m sufficiently large, each of the ratios 



, , • • • • 

u u ^ u , 



may be made to diflfer from I as little as we please : and there- 
fore also the geometric mean between these ratios, viz. 









T^^Jf l>y sufficiently increasing n, be made to differ from / by 
a quantity B as small as we please ; we shall have, therefore, 

1 1 

Consequently, supposing n infinite, in which case S vanishes, 

limit of (uj* = / ; 

1 
but limit of (ujT = L, 



^ XXX x*" 

Ex. 1. 6*=!+-+ — + + ....+ +. . . 

1 1.2 1.2.3 1.2...W 



-^^ = ultimately when n is infinite ; 



u « + I 



n 



the series, therefore, is convergent for all values of x. 



CONVERGENCY AND DIVERGENCY. 79 
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«n.l 
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x^ 


2)' 


= 0, 










u 

n 


(2w + 


1) (2w + 





when n is sufficiently large. Therefore the series for cos x is 
always convergent. Similarly, the series for sin x is always 
convergent. 

Ex. 3. The series for 

log(l + a;) = :i;-- +^ +(-l)«-»?!+ 

^ ^ ^ 2 3 ^ ^ n 

is convergent or divergent, according as ar is less or greater than 
unity. The same is true of the series for tan'* x^ viz. 

^ x' 
^" 3"'6 ■•••• 



Ex. 4. (l^xf^l^pLX-^: (iOi_JL) ^ 

. fi(fx- 1) (i^-y^4 1) ^, . 

1.2.. ..w 

/. _2iJ = ^ X =: - X ultimately when n is infinite ; 

u^ « + 1 

and therefore if a: is less than 1, the series is convergent. 



93. (3) The series w, - w, + Wg ~ w^ + &c. 

is convergent, if the numerical value of its terms decreases 
without limit. For we may write it in either of the forms 

or (Wi - u^) + K - W4) + ("4 - ^5) + &c-> 
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which shew that the sum of the series is greater than u^ - u^, 
and less than u^ : hence it is convergent. 

94. There is another test which it may be as well to notice ; 
it is as follows : the series 

w, 4 w, • • • • + ^n "^ • • • • 

is convergent, or will become so, if the inferior limit of v f\ 

n 

be greater than 1, when n is infinite; and divergent in the 
contrary case. The proof is very much the same as that of the 
test in Art. 91. 



96. If the general term be imaginary, and therefore of the 

form 

w = c + w?^ V- 1 , 

ft n n ' 

we must assume 

v^ = Pn cos 6^ and w^ = p^ sin »„, 

Then p„ is greater than t?^, and also greater than w^. Also 

K = Pn (cos e, + V- 1 sin flj, 
and we may consider the series to be made up of two, viz. 

(«?, + t?j + .... H- e?^ +....) + >/- 1 (t^j + t<?2 + .... + «^^ + ... .)> 

and the tests may be applied to them separately. Both will be 
convergent if the series of moduli, 

Pvp2 Pn^ 

form a convergent series. 
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On Interpolation. 

96. The object of Interpolation is to insert amongst the 
terms of a given series new terms subject to the same law as 
those ofthe given series. 

The practical use of this part of the subject is seen in cases 
where a series of values of a quantity being obtained, either 
from calculation or observation, it is necessary to insert between 
them other values, such as would have been obtained from 
similar observation or calculation without actually performing 
these. 

The function which expresses the general term of the series 
is therefore what is required to be known: and it is to be 
found from a limited number of its numerical values. This 
caniidt, however, be effected from these data, but if the terms 
given are sufficiently numerous, and differ but little from each 
other, the function can be found to a great degree of accuracy. 

97. We proceed to consider the most usual cases of inter- 
polation. 

(1) Having given n values of a function y*(a;), 

^i> ^a> ^8> • • • • ^«> 

corresponding to n values of the variable 

o, a + A, a 4- 2A, . . . . o + (n - 1) A ; 

to find any intermediate value, Bjs/(a + A). 

Let k = mh. 

Then, •/ 

f(x + mh) = f(x) + mAf(x) + ^?^^) ^y(^) + . . . . 

k 
if we make x = a, and put r- for m, since the value of /(x) 

in this case is u^j we have 

G 
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But AWj = ^2 - ^P 

&c. = &c., 
j^n-i f ,N (n-l)(w-2) 

^ X = ^'n - (» - ^n-l + Y72 ^"-^ " •" * "'' 

Hence these quantities Aw,, A'w,, .... are known, and on sub- 
stituting them in the series for f{a + A), we have the value 
required. 

98. In using this formula, the easiest way of finding the 
values of Aw,, A'w,, &c. is not from the above expansions, but as 
follows: write down the series of values Wj, ^, «g. . . .w^ one 
under the other, and subtract each from the one under it, 
arranging the results one under the other in a separate column ; 
the first of these is AWj. Perform the same process of sub- 
traction with the terms in this column, and write the results 
as before, the first of these is A'w,, continuing the operation 
as long as there are any terms. 



Ex. Having given the values of log,ol23, log,ol24, logjol25, 
log,ol26; to find log,o 123.25. 

If a = 123, A = .25 = J, and u^ - logj^a:, 

w. = 2.0899051 

t/^^, = 2.0934217 



u^^^ 2.0969100 
w„^3= 2.1003705 



A 




35166 


A' 


34883 


- 283 


34605 


-278 



5. 



Now 



7 A *(*-!) A2 A(A-1)(A 



-2) 



1.2 



1.2.3 



AV 



terms after these being neglected. 
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Putting for k its value in the second member, 

.-. log (123 + J) = 2.0899051 

+ .0008791 
+ .0000027 
= 2.0907869. 

99. (2) Having given a series of consecutive and equi- 
distant values of a function, where one is deficient to insert 
that one. 

Let Wp w,. . • •«„ be the values of the function, corresponding 
to the values 

o, a + A, a + 2h, . . .a + (w - 1) A 

of the variable. In this case the terms given must be so 
numerous, and near each other, that A**'X is zero, or nearly so : 
this will almost always be the case in tabulated results. 

^ ,N (7i-l)(w-2) 
••• ^" «^, = «^n - (w - 1) K^i + ^^ f^ K-2 +. . . . ± Wi = : 

from this equation any one of the values, as w^, may be deter- 
mined. 

If two values are deficient, we must make another sup- 
position, viz. that A""'«i = 0, and we have thus two simple 
equations for the determination of the two quantities. 

Ex. Given the logarithms of 121, 122, 124, 125; to find 
that of 123. 

If a = 121 and u^ =* logio^> 

W« = logiol21 = 2.0827854, 

u^^ = logiol22 = 2.0863598, 

w„,3 = log,J24 ==-2.093421 7, 

u„^^ = logiol25 = 2.096^100; 
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.-. putting A'Ua = u,^ - 4w^, + 6w«,, - 4u^^^ + t4. = 0, 

= J {4 (4.1797816) - (2.1796954)} 
= J {16.7191260 - 4.1796964} 
= J {12.5394306} 
= 2.0899051. 

100. To interpolate any number of equidistant terms between 
every two consecutive terms of a given series. 

Let the given series be w^, w,. . . .w^, and let (r - 1) equi- 
distant terms be required to be inserted between every two 
of them. 

The series will then become 

r r r r r 

Now u.^P = w, + ^ Au, + ^' A"o ■ A\ + 

and therefore by making /> = 1, 2, ... .(r - 1), in succession, 
since the values of Aw^ A'Wj. ... are known from the given 
series, we get (r - 1) terms between u^ and u^. Next, taking p 
from r + 1 to 2r - 1, we get (r - 1) terms between u^ and w„ 
and so on. 



Lagrange's Formvia of Interpolatum, 

101. By the aid of this formula the function can be found 
which for any number (n) of values a^ a^ . . . . a^ of the 
variable, proceeding according to no regular law, assumes 
the values 

Wj, WgJ ^s> • • • • Wn» 



INTERPOLATION OF 8BRIBS. 85 

Let w, = -4j + wdj + . . . . + -4„, 

where A^y A^y .... A^ are each functions of a: ; the form 
of u^ is in no wise restricted by this assumption. 

Then, when x =» Oj, the second member must reduce itself 
to t^ J, in order to which all its terms must contain the factor 
X - a,, except one, as A^ for instance. Similarly all the 
terms, except one, as A^ for instance, must contain the factor 
X - a^y and so on. Hence each of the terms will contain all 
except one of the factors a; - a^ a; - a^, . . . . a? - «,, and that 
one different for each term. 

Whence we gather that the form of w. is, 

u, = Al(x - a;){x - a^...{x - a^ + A^{x - a,){x - a;)...{x - aj, 
"t" ■ • • • 

But when x = a^y u^ becomes Wj, and hence 
Similarly, 

&c. = &c., 

^ (a; - aj(a? - g,) (a: - gj] ^ 

(«a - «i)K - ^a) • • • • («2 - O '* 

+ &c. 

, {x - a^{x -a,) (x- a^,,) ^^ 



86 INTERPOLATION OF SERIES. 

Ex. To find the function of x which^ when x - ^, Zy 2, I 
shall assume the values, 5, 6, 7^ 8, 

^ _^ (x- 8)(a: - 2)(x -I) ^ (x- 4)(a; - 2)(a; - 1) 

1.2.3 1.2 

(x-4)ix-'S)(X'l) __ 3 (a;-4)(a?-8)(a;-2) 
2.1.1 3,2.1 
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1. A 



EXAMPLES. 
I. 

Direct Method of Difference*. 



8. A2.^ = 2.^. 
of - I or + I 

8. A2-8in|:=2-8inA^sinljJ. 

4. A(2-rinlJ = (2--8mAJ. 

5. Atan- 5-. 

6. A2-tani = -2-tan|j^tanliJ. 

7. A (- 2)- sin I = (- ir a-" sin A (cos A J, 



8. A— L- = -2 



. e 

em — r 

2" sm — 2' em — 

2" 2' 

tan — ; 
1 "^2"' 



9- A g- = -2ST 

2*tan -- 
2* 
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tan Ay 
10. A/ ' jiV"-^-V %~)- 



' ^ \2' tan l) 






cos- sm^ 



11. A J All " . ^ , !, 



F4J"( 



2**'cos-^j 



2«+i i 



12. Acot2-e = --T ^ 



•♦i/i* 



sin 2**'e 
jg log (2 sin 2^0) ^ ^ log tan 2^0 

j4 ^ 1 1 

(2- Bin -J (2-^co8^" 

15. A2' tan-^ - = 2' tan'^ f ^ , ^^ ^^J ) . 

2* \4.2'" + 3.2'6'/ 

16. From the expression 

A* sin «6 = f 2 sin — j sin {xO + - (?r + hffj}, 

obtained on the hypothesis of an indeterminate increment h, 
deduce the value of — -z . 



Differences of Zero. 

1. Investigate the following relations : 

A'O* - A'-'O' +...± AO' T 0* = 1 ; 

AO- A'O- A'O* ^ 

and r- ±...± = 

12 X 
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2. Shew that 

{log(l + ^)yf(A)0' = x(x - !)...(« - n + 1)/(A)0— ; 

and that 

/{(I + A)*} 0- = n-/(l + A) 0-. 

S. Prove that 

{/(l+A)+/(^)}0--«-0; 

and |/(1 + A) -/(]^)) 0" = ; 

and thence deduce the expression 

AO*' A'O*" A'O'" A**©** 

""2 ¥' ■*■ "F* "^'"^ 'W " ^* 

4. Assuming /(I + A) in (3) to be 

J^^log(l^.A). 

prove that 

AO** A'O** A'"0'* 

in all cases^ except where z = I. 

5. If /j(A) and /3(A) denote any two functions of A, then 
the following equation will hold good^ 

{/(A)/,(A)}0- -/,(A)/,(A'){0 + 0'}-; 

where in the second member the unaccented A is to be re- 
ferred to the unaccented powers of 0, and the accented A to 
the accented powers. 
Hence shew that 

{(1 + Ar/(A)} 0- -/(A)(n + oy. 

6. Prove the following formula for changing the inde- 
pendent variable from a: to e*. 



( 



dY AO" d AV , (f A'O" , d' 
^ dy) ^~r^ dy'^ 1.2 ^ df'^ 1.2.3 ^ df '^'" 
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7, Shew that iSu^^ a + hx, 

8. Investigate the expression^ 

1.2...(« + ») 1 A J 
and interpret its meaning. 



II. 
Inverse Method of Differences. 

1. Find S i-jT^ ' ^^ 2(** + aj* + 1). 

2. Find 2 Z^-. — "^ , and 2 ( ? \ . 

8. Find 2a"sin':t;9^ and shew that a particular value of 
2*(a""co8a?6) is 

oT' (P-^ ' ay cos xO 
(1 - 2a cos e^-ay 

4. Investigate 2 -7 — —7 — -^ — -r and 2 t--^ — rr^- r; • 

5. Find 

2 (2a: + 3)', 2(«" + x){2ix? - 6/, and 2(a;* + 1)». 

6. Reduce u^ju^ to a series of terms^ each of which is the 
product of consecutive terms of the arithmetical progression of 
which Ug^ a -¥ hxia the general term. 

7. Apply the same process to (w,)'(w.^,/. 
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Calculus of Oenerating Functions, 
1. Shew that 

1 - 2^ cos 6 + ^ 
and deduce O singed. 



2. Prove that 



^/*rT/-s^?/*a- 



2 log 

ft 

3. Determine the generating functions of log z, a^, and or**. 

4. Shew that 

Qe"^Jdxd^^ » (a - \ogt)'^Ou,, 

and Gc^-^^a-'u^ = (1 - «^)** ^<?w.. 

5. Find t<^ from the following relations : 

1 - V(l - 40 

6. Having given that Ou^ « F(t), shew that 

6? (r + 2'» +...+ aj») w. = ii^ {JF(^ + ^A) - F(t)} Q\ 



Numbers of BemouiUi. 
1. Prove that B^,^ being the z^ of these numbers. 



2 + A 2« 

^. Also shew that 
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III. 

EqtmUoiu of Differences. 

1. Solve the following equations : 

tt^i - (« + l)«. « 1.2. ..(a: + 1); 
and (x + l)"(w.n - «wj « a*. 

2. Integrate tt^,i - cw, « sin 3:6, 
and w^i - ca?w, « cosa;0. 

8. Solve 

<*«.« - «(« + 2)tt,,, + b(x + !)(« + 2)w. - 0; 
and tt^ + ap^u^i + ^^'w, = ^^P'^' 

4. Find the integral of 

tt^, + a0 (a; + l)tt^j + J0(a;)0(2? + l)w. - c, 
and that of Au« + A't«« « 42u^. 

5. The solution of 

/ oM n (« + 2)' (x + l)»(a: +2) 



IS 



• * « + 1 * v^+ 1/ ' 



6. If (0"-4(uJ{(f/J'-H}, 

7. Solve (wj - 4uu^, + a(2w^ + u^,) = ; 
and (uj + 2w.w^i - a(4w. - u^,) = 0. 

8. Prove that y = oa: + a', 

and y = U-»(-in'-7. 

4 

are distinct solutions of the same equation in finite differences. 
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9. Find the indirect integral of 



Au. 



10. Solve u^ - 6w^j + 1 lu^i - 6fi^ = 0* + «* ; 
and u^i - x^u^ = (1.2. ..a?)*. 

11. Shew that the general integral of the equation 
is comprised in that of the equation 

and is u^ = A + Ba' + Ca'" ; 

a being one of the imaginary cube roots of unity^ and A, B, C 
being connected by the equation 

A'-^ff+C^SAia + BC). 

Also shew that 

u =3 va tan C, cos + CL sin ) 

• V ' 3*3/ 

satisfies the equation. 

12. If w..iW» - «.(w«.i - t/J + 1 = 0, 

w = tan ( (7 + S tan"* — j . 



Simtdtaneous Equations. 
1. Find w, and t?^ from the following equations : 

tt^, - 2t>. = a«. 



^^i + 2w. = a-. 



2. Solve aw.^, - J (t/^ + t>J = 0, 
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8. Solve 


«^. - 2o„, + «. = 0, 




».« - 2«^. + t». = 0. 



4. Solve tt^ - 8t>. = 27w.,i, 

5. Solve w^-C^'^i+O"*^.* 



Equations of Partial and Mixed Differences. 

1. Find u^^^ in each of the following equations : 

2. The solution of the equation^ 

«^, = i^.^...^ [i>r^. y3^^ + <^»-) ' 

where P^A^^ denotes the product of all the successive values 
oi A^^ obtained by giving x all successive values up to x^ and 
C^^ is an arbitrary function of y + a:. 

3. Integrate the equation 



IS 



«-.»-^«-r=0- 



4. The solution of 



d 
A u ^a -r- u • 

• ••* dff '*^ 



is ^... = ^'^'"(^J*(y)' 

denoting an arbitrary function. 
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where a and /3 are the roots of «* - a« + J = 0. 
6. The integral of 

A » A ^ T ^ 



IS u ^ a' e 



where a and /3 are the same as in the preceding. 



Functional Eqtcations. 

1. Find a particular solution of the following : 

2. If /(»)+«/(-«) = «", 
8. If /(a:) - o/^ L) = «-, 



e* + ae' 



•^(*) = Tiv 
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4. If (/(,))../^l^) . A. 

/« - (HI "-^ 

5. If /(I + X) +/(l -«)=!-«', 



IV. 
On Seriet. 

1. If 1.2.4 + 2.8.5 + 3.4.6, to x terms, = S^ 

, „ ar (« + 1) (a: + 2) (8« + 1 a) 
then A, = — 

2. If 1* + 8' + 5*... + (2* - 1)* = S,, 

a;(4a:'-l)(12a^-7) 
then o^ = — — . . 

1 1 



S. If 



V2(l + V2) (1 +V2)(2 + V2)"(2 4-V2)(3 + V2) 

+ • • • • to .6 termSj «= ^^^ 

then ^. 



• V2 (« + V2) ' 



4. If 1.3' + 3.5' + 5.7' + ... to X terms, = S^y 

« a:(6a;'+ 16aj' + 9a;- 4) 
then 5. = -^^ ' . 



5. If + — z — z — :: — + — to x terms, - S ^ 

5.7 7.29 29.79 79.245 ' -' 

then ^^-TTz^ 



• 20 4 {2 + (- 3)-*} ' 
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6, The series whose (x + IJ^ term is 

(p + qx) ^ 



^.'^^i-^^m-i 



, where «^ = a + Aa:, 
has for its sum 

S^ g ( ^ l—\; 

when the following equation of condition holds good^ 

p _a m-l 
q'^h't'l' 

l' 4 2'.4' 3'. 4' 

7. Thus if — ^ + -^— + — ^ to a: terms, = S. 

2.3 3.4 4.5 

*v o 2 4 a: - 1 ^. 
then /S = - + , 4*. 

• 3 3 a:+ 2 

8. If 1 cos' 6+2 cos' 20 + 3 cos' 36 +... + « cos'a:6 = iS , 

^, cv a:.(a: + 1) a: sin (2x +1)0 1 /sin xOY 

then iS. = —4 -' + - — ^ . n - 7 -^-s" • 

4 4 sinO 4\ sin fly 

9. If sin 26 COS0 + sin46 cos 30 +...+ sin2a:6 cos(2a:- 1) ^ = S^^ 

. ^ sina:(0+0) . .^ .^ .^ 8ina:(6-0) . ^.^ .^ ., 
then^y^^: , . }^ ZU m{e-{'X(e'¥4>)]-^ ^ . ^ ^; sm{6+a;(6-0)}. 
2sm(6+0) *• ^ ^'^' 2sin(6-0) *■ ^ ^'^^ 

10. If 1' sec' 6 + 2' sec' 20 + 4' sec' 46 + ... + to ar terms, = S,, 
then S^ = 2" cosec' 2^6 - cosec' 6. 

11. Apply the expression in Art. 23 to sum the following 

series : 

1* - 2" + 3** - 4* + to 00 , 

1* + 2" + 3* + 4*» + to 00 , 

,, 2'»A 3"A' 4*A' 
1 1.2 1.2.3 

1« « 3« + 5- „ 7» + to 00 . 
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Hence shew that 

1-1 + 1-1. ...tpooB J. 



12. If ^^j «_-_+_-.... to 00, 

(- M-T 



Becurring Series. 

1. Smn to X terms the following recurring series : 

1 + 13 + 49 + 157 + .... 
whose scale of relation is 1 + /> + ;. 

2. Sum also the following : 

2 - a - a* + 2a' - a* - a' + 2a" - a' - a' + ... . 

S. The numbers A*0% A"0**', A-0***... 
form a recurring series whose scale of relation is 

-4j - -4j + -4, - ... ± ^^, 
where -4, = 1 + 2 + 3 + ... + », 

-4, = 1.2 + 2.3 + 1.3 + ... +(»-!) », 



-4_= 1.2.3 ... n. 



Interpolation. 
1. Having given 

Z tan 46° = 0, i tan 45°.30' = 0.0075803, 

Z tan 46° B 0.0151628, L tan 46°.30' = 0.0227500; 

find L tan 45°.20'.20". 
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2. Given log^^ 510 = 2.70767018, 

logio511 = 2.70842090, 
log,o513 = 2.71011737, 
logio514 = 2.71096312; 
it is required to insert the deficient value 

log,o512 = 2.70926996. 

3. Given t?o, Uj, t?^, v^; required t?, and t?,, 



4. Given 



Vj- - 




10 


Mm — 


2t?o + 


5v^ + lOo^- 3t?g 


t?^« - 




10 


tan- 


^10« 


1.471127674, 


tan" 


Ml = 


1.480136439, 


tan" 


^13 = 


1.494024435, 


tan" 


^15 = 


1.504228163; 


tan"' 


11.63 


= 1.485022707. 



shew that 

5. Find the function which, when x assumes the values 
2, 4, 6, 8, shall assume the values 9, 10, 11, 12. 



C 100 ) 

Table of tjie numerical values of i"0" for all values of 
m aud n &om 1 to 10. 
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